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The fundamental assumption of statistical mechanics is that the system is equally likely in any of 
the accessible microstates. Based on this assumption, the Boltzmann distribution is derived and the 
full theory of statistical thermodynamics can be built. In this paper, we show that the Boltzmann 
distribution in general can not describe the steady state of open system. Based on the effective 
Hamiltonian approach, we calculate the specific heat, the free energy and the entropy for an open 
system in steady states. Examples are illustrated and discussed. 
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INTRODUCTION 

The fundamental postulate in statistical mechanics is 
as follows: For a conservative system the path of the 
system in phase space passes through all points of the 
energy surface and in such a manner that this surface 
is covered uniformly. This postulate (called ergodic hy- 
pothesis) indicates that a conservative system in equilib- 
rium does not have any preference for any of its available 
microstates, i.e., given n microstates at a particular en- 
ergy, the probability of finding the system in a particular 
microstate is p = By using the ergodic hypothe- 

sis, one can conclude that for a system at equilibrium, 
the thermodynamic state which results from the largest 
number of microstates is the most probable macrostate of 
the system. With these results, the probability pi that a 
macroscopic system in thermal equilibrium with its envi- 
ronment in a given microstate with energy Ei can be de- 
rived. Pi = exp{- PEi) / e.xp{- PEj) and P = l/ksT, 
which is the so-called Boltzmann distribution. 

It is believed that a generic system that interacts with 
a generic environment evolves into an equilibrium de- 
scribed by the Boltzmann distribution Experience 
shows that this is true but a detailed understanding of 
this process, which is crucial for a rigorous justification of 
statistical physics and thermodynamics, is lacking^-Q- 
The key question is to what extent the evolution to equi- 
librium depends on the details of the system-environment 
coupling [8|. In fact, detailed analysis shows that it is 
not always the case that the system evolves into the 
equilibrium state described by Boltzmann distribution. 
Then two questions arise: (1) Given a dynamical process, 
how far does the steady state differ from the equilibrium 
state? (2)What is the free energy and specific heat of 
such a system? 

Here, we shall answer these questions by considering an 
open system coupled to two independent environments at 
different temperatures. The dynamics of the open system 
is assume to fulfill a master equation in the Lindblad 
form, steady state is achieved by the effective Hamilto- 
nian approach. We find that the steady state is not the 
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FIG. 1: (Color online) A given system immersed in two envi- 
ronments at temperature Ti and T2 , respectively. Eventually, 
the system arrives at a steady state depending on the system- 
environment couplings and the temperature. 



thermal equilibrium state even if the environments have 
the same temperature. Free energy, the specific heat and 
the entropy are calculated and discussed. 



FORMALISM 

Imagine to prepare an open quantum system sur- 
rounded by two independent environments at different 
temperatures Ti and T2, respectively. See Fig IT] In the 
weak coupling limit and under the Markov approxima- 
tion, the dynamics of the quantum system is govern by. 
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[H, p] + Cip + C2P, 



(1) 



where Cip = ^J2k^kiTi){2LkpLl - LlLkP - 
pLiLk),C2P = ^j:MT2){2XkpXl - xlXkP - 
pXj,Xk). p{t) represents the reduced density matrix of 
the system, H represents the free Hamiltonian of the sys- 
tem and Lfe(ifc) and Xk{Xl) are operators of the system 
through which the system and its environments coupled. 
The master equation can be solved by the effective Hamil- 
tonian approach 0. 

The main idea of The effective Hamiltonian approach 
can be outlined as follows. By introducing an ancilla, 
which has the same dimension of Hilbert space as the 
system, we can map the system density matrix p(t) 
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to a wave function of the composite system (system 
+ ancilla). A Schrodinger -like equation can be de- 
rived from the master equation. The solution of the 
master equation can be obtained by mapping the so- 
lution of the Schrodinger -like equation back to the 
density matrix. Assume the dimension of the Hilbert 
space for both the system and the ancilla is N, and 
let |£^„(0)) and |em(0)) denote the eigenstates for the 
system and the ancilla, respectively. The mathematical 
representation of the above idea can be formulated as 
follows. A wave function for the composite system in 
the A/'^-dimensional Hilbert space may be constructed as 

p{t) ^ \^{t))p = j:l,n=lPmn{t)\E^iO))Km, whcrC 
Pmn{t) = {EUOMt)\Enm. NotC that p{^\^) p = 

Tr(p2) < 1, i.e. this pure bipartite state is not normalized 
except when the state of the open system is pure. With 
these definitions, the master equation {H = 1 hereafter) 
can be rewritten in a Schrodinger-like equation^ 

«^l*W>P=^cff|*W)p, (2) 

where Hcff is the so-called effective Hamiltonian and is 
defined by 

T^eff — ^ — 

+^J2^k{T^)L^Lk + ^Y,^k{T2)X^X,, (3) 

fc fc 

where H = H - ^Y^k IkiWlL, - f Y^k rkiT2)xlXk, 
and X^ are operators for the ancilla defined by 
(O = L,X,H), (e„(0)|O^|e„(0)) = (£;„(0)|Ot|S™(0)). 
By this effective Hamiltonian approach, it is easy to prove 
that the steady state ps can be given by mapping the 
eigenstate of T-LcS with zero eigenvalue. Namely, calcu- 
lating |i?o) by 

H^s\Ro) = 0, (4) 

we can obtain elements of the steady state density ma- 
trix, ps.mn = {E^mps\Enm = (i^™ (0) | (e„ (0) | i?o> • 
Given a steady state, the single-particle energy U would 
be equal to, 

U = Ty{psH). (5) 

If a system consists of many non-interacting particles, 
the total energy equals to the sum of the single-particle 
energy. The specific heat for a single particle now would 
be given by 

Given the steady state density matrix p5, von Neumann 
defined the entropy as 

<5' = ~P5lnps, (7) 




FIG. 2: (Color online)The specific heat, the free energy and 
the entropy for the open two-level system as a function of 
temperature (in units oi The parameters chosen are 

7 = 0.2S1, r = 0.3f2. The energy was plotted in units of 
hO,. The units for the specific heat and entropy were set 
accordingly. 

which is a proper extension of the Gibbs entropy to the 
quantum case. We note that the entropy S times the 
Boltzmann constant ks {ks = 1 in this Letter) equals the 
thermodynamical entropy. If the system is finite dimen- 
sional, the entropy describes the distance of the steady 
state from a pure state. 



EXAMPLES 

To illustrate the general formalism, we present here 
three examples. In the first example, we consider a 
two-level system coupled to two independent environ- 
ments at different temperatures Ti and T2, respectively. 
The dynamics is described by Eq.(IT]) with 71 (Ti) = 
7ni, 72(7"!) = 7("-i + 1)> Li = a+, L2 = cr^ , 
ri(T2) = r(2n2 + l), Xi = and n, = ^^^{i = 1, 2). 
The system Hamiltonian is specified to be i7 = ^cr^. By 
the effective Hamiltonian approach, we arrive at density 
matrix elements of the steady state ^ = ^^'1'!^^ with the 

Poo T2|±,i 1 

trace preserving condition pn+poo = 1- Figl2]depicts the 
specific heat Cti , j the free energy U and the entropy 
5 as a function of temperature Ti and T2 . We note that 
for Ti — > and Ti — > 00 the specific heat Cti tends to 
zero, the population is then said to be frozen. For a fixed 
Ti being of order of Q, Cti decreases as T2 increases, and 
Cti tends to zero as r2 — > 00. behaves similarly. The 
free energy U and the entropy S approach constants with 
Ti and T2 tend to infinity, confirming that the population 
is frozen at sufficiently high temperatures. At low tem- 
perature, U and S increase as the temperature increases, 
indicating that the degree of mixture of the steady state 
grows with the increasing of temperature. 

We take two coupled qubits subject to decoherence 
as the second example. The system Hamiltonian is. 
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FIG. 3: (Color online) The specific heats Cti and Ct2 versus 
temperature Ti and T2 for two interacting qubits dissipatively 
coupled to two independent environments. f2i — , 7 ~ 
0.2r2i, J — 0.2r2i,r = O.Sfii are chosen for this plot. 
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FIG. 4: (Color online) Cti and Ct2 as a function of tem- 
perature and the coupling constant J. T2 = 1.5 for (a) and 
Ti = 1.5 for (b). The other parameters are the same as in 
FigEl 

H = Sli|e)i(e| -I- r22|e)2(e| + Jcrfcr2- Suppose that the 
qubit 1 interacts with its environment at temperature Ti 
via (Tj~ and , while the qubit 2 through cr| couples 
to its environment at temperature T2. The Liouvillian 
superoperators are then, 

+ -'yni{2a+pa^ - pa^at - crj;(7+p) (8) 

and 

C2P - ^r(2n2 + l)(2af paf - 2p), (9) 

where ni = ^371^731; ^2 = pj^TT^- ^- as Ti and 
T2 tend to 00 as FiglHJb) shows. For Cti , however, it 
approaches a constant as T2 00 when Ti takes a value 
of order of fi^ (FiglSI^a)). At Ti = 0, Cti is always zero. 
For fixed Ti, Cti tends to a constant with J — > 00, while 
always tends to zero as J ^ 00, as FigHlshows. Equi- 
librium statistical mechanics tells us that the population 
of excited states (given by the Boltzmann distribution) 
grows as the temperature increases. The populations ob- 
tained from the steady state are different (see Fig[5]). E. 
g., the population of the ground state increases as the 
temperature T2 increases (FiglS] (left)), whereas the pop- 
ulation of the second excited state (labeled by (3)) de- 
creases as the temperature increases. Similar observation 
can be found from FigOl^right), the population of the first 



FIG. 5: Probability for finding the system in the eigenstates of 
the Hamiltonian. Qi = ^2,7 = 0.2^1,7 = 0.2^1, F = O.Sfii. 
(l)-(4) label the eigenstates by the corresponding eigenvalues 
in increasing order, i.e., (1) labels the lowest eigenstate, while 
(4) the highest eigenstate. 

excited state decreases as the temperature Ti increases. 
We will quantify this difference between the steady state 
and the equilibrium state by fidelity later. 

In the third example, we consider a damped har- 
monic oscillator. The system Hamiltonian takes H = 
wa^a. Consider a simple system-bath (at temperature 

Ti) Hamiltonian of the form Hint — J2j dji^j'^^ + ^]'^)- 
The damping rates follows by the standard procedure, 
liiTi) = 7^1, Li = at, 72(Ti) = 7(711 -I- 1), L2 = a. 
Suppose that the harmonic oscillator interacts with the 
environment at temperature T2 through (a^-l-a), the cor- 
responding damping rate is ri(T2) = r(2n2 + 1), and 
Xi = {a^ + a). All these together give a master equation 
for the damped harmonic oscillator, 

d 1 
ig^P [H, p] + -ji{2a'' pa - paa'i - aa'' p) 

+—j2{2apa^ — pa^a — a^ap) (10) 

+ ^ri[2{a'^+a)p{a''+a)-{a''+ayp-p{a''+a)'^]. 

Figini shows the specific heat, the free energy, and the 
entropy as a function of temperature. We note that the 
specific heat Cti {CT2) is vanishingly small for Ti — >■ 
0{T2 — > 0), it rises rapidly when Ti (T2) is of order of 
uj and approaches a limiting value, which depends on 
the damping rates Fi and 7. Note that at sufficiently 
high temperature, Ct^ {i = 1,2) are the same as that in 
equilibrium statistical mechanics. 

The results presented in the examples clearly show that 
{i — 1,2), U and S behaves different from those given 
by equilibrium statistical mechanics. The differences re- 
sult from the deviation of the steady state density matrix 
from the equilibrium thermal state (Boltzmann distribu- 
tion). We will use the fidelity to quantify this deviation. 
Fidelity as a measure of distance between two states is 
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FIG. 6: (Color online) Illustration of Cti , Ct^ , U and S as a 
function of temperature Ti and T2. In this plot, the dimension 
cutoff is = 100, and 7 = O.Soj, F = 0.2a;. The energy is in 
units of huj and the temperature is in units of 




with Fi = can describe the thermalization of the oscil- 
lator, the simulation shows that this is exactly the case 
(FiglTjc)). Differently, for the coupled two qubits system, 
the steady state given by the master equation (Eqs.® 
and ^) with F = is not the equilibrium thermal state 
at low temperature (see Fig 13 (d)). 

One may have doubts about the realizability of the 
master equations Eqs. and dTO]). The tech- 

nology in engineered reservoirs j 13. 14] shows that there is 
no problem to simulate an reservoir in which the system- 
environment coupling and state of the environment are 
controllable. 

In summary, we study the deviation of steady state 
from equilibrium thermal state for an open system. The 
specific heat, free energy, and entropy are calculated and 
discussed. This study applies to several occasions, where 
a great many of physical phenomena of interests concern 
collective behavior of an open system in steady state. 
This work provides the exact solution for the nonequi- 
librium distribution and statistical quantities for steady 
states, thus giving insight on how to build a statistical 
mechanics for open systems. 

This work is supported by NSF of China under grant 
Nos 10775023 and 10935010. 
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FIG. 7: (Color online) Distance between the steady states 
and thermal states (Boltzmann distribution) measured by the 
fidelity as a function of temperature, (a) and (c) are for 
the damped oscillator, while (b) and (d) for the two coupled 
qubits. The steady state in (a) and (b) are obtained by solv- 
ing the master equations Eqs. ((8|, ((9]) and (|10p . The steady 
states in (c) and (d) are for Fi = and F = 0, respectively. 
^1 = ^2, J = 0.2^1. 



an important concept in quantum information theorv[lo|. 
The well-known quantum fidelity for two general mixed 
states pi and p2 is given by the Uhlmann's fidelity 11 1 



F{pi,P2) = Tr J y/plp2y/pi, 



(11) 



this fidelity possesses many advantages such as concavity 
and multiplicativity under tensor product and it satisfies 
all Josza's four axioms fl2j. FiglT] shows the fidelity be- 
tween thermal states and the steady states. For both the 
coupled qubits and the damped harmonic oscillator, the 
fidelity arrives at its maximum when the temperatures 
tend to infinity (see FiglTja) and (b)). Note that the 
steady states depends on how the system couples to its 
environments. For example, the master equation Eq. ljlOp 
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